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Accurate spacecraft tracking by the NASA DSN antennas must be assured during 
changing weather conditions. Wind disturbances are the main source of tracking 
errors. This article presents the development of a wind- force model and simulations 
of wind-induced pointing errors of DSN antennas. The antenna model includes the 
antenna structure , the elevation and azimuth servos , and the tracking controller. 
Simulation results show that pointing errors due to wind gusts are of the same order 
as errors due to static wind pressure and that these errors (similar to those of static 
wind pressure) satisfy the velocity quadratic law. The methodology presented is 
used for wind-disturbance estimation and for design of an antenna controller with 
wind-disturbance rejection properties. 

I. Introduction 

Reliable estimates of the pointing error due to wind 
disturbances are required for operating antennas as well 
as for designing effective controllers. There are two types 
of wind-induced pointing errors: steady-state error due to 
static wind pressure and dynamic error due to wind tur- 
bulence. Steady-state error models of antennas have been 


1 Independent consultant to the Ground Antennas and Facilities En- 

gineering Section. 


investigated by Katow and McGinness [1,2]; the dynamic 
error of the DSS-14 antenna has been analyzed by Mas- 
soudi [3]; and the 100-m NRAO antenna has been analyzed 
by Hill. 2 


2 R. Hill, “Servo Design Studies and Pointing Performance Analysis 
for a 100-m Aperture Radio telescope,” in Servo Design Specifica- 
tion for the 100-meter Unblocked Aperture Telescope for the Na- 
tional Radio Astronomy Observatory , JPL D-7477 (internal docu- 
ment), Jet Propulsion Laboratory, Pasadena, California, Septem- 
ber 1990. 
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In this article, a recently developed antenna model [4] 
has been used to analyze wind disturbances. The struc- 
tural finite-element model, as a part of the antenna model, 
includes modes of up to 10 Hz. The dynamic wind pres- 
sure consists of the spatial pressure on the antenna struc- 
ture and the wind-force time history. In the new approach 
presented in this article, the wind force is modeled as a 
pressure on the antenna dish rather than as an equiva- 
lent elevation torque commonly used as a replacement for 
wind action. Wind-pressure distribution is obtained from 
the J PL/IDEAS finite-element code, 3 4 5 based on wind tun- 
nel measurements [5]. 4,5 The wind-force time history is 
another new addition. It is a linear filter that shapes the 
wind-force time profile obtained from the Davenport wind- 
velocity spectrum [6]. 

The results obtained consist of estimates of pointing 
error due to static wind pressure and wind gusts for 60- 
and 90-deg antenna-elevation angles. Results show that 
the pointing error due to wind turbulence is of the same 
order as the error due to static pressure, thus it must not 
be neglected. 

II. Wind-Velocity and Wind-Force Models 

The wind velocity v is a combination of a steady-state 
flow (mean velocity, or a quasi-steady component, see [10]) 
v m , and a turbulence (gust) v t 

V = v m +v t ( 1 ) 

The component is a random process with zero mean 
and standard deviation cr v . The standard deviation and 
other wind parameters are determined from the Davenport 
model [5,7]. In this model, the wind- velocity spectrum S v 
for the wind at a height of 10 m is as follows: 

S.(/) = 4.?/- , fT ^i7J (2) 


3 R. Levy and D. Strain, JPL- IDEAS Finite Element Analysis and 
Design Optimization Program, JPL NPO-17783 (internal docu- 
ment), Jet Propulsion Laboratory, Pasadena, California, October 
1988. 

4 R. B. Blaylock, “Aerodynamic Coefficients for a Model of a Para- 
boloidal Reflector Directional Antenna Proposed for a JPL Ad- 
vanced Antenna System , ** JPL Memorandum CP-6 (internal doc- 
ument), Jet Propulsion Laboratory, Pasadena, California, May 1, 
1964. 

5 N. L. Fox, “Load Distributions on the Surface of Paraboloidal Re- 
flector Antennas,” JPL Memorandum CP-4 (internal document), 
Jet Propulsion Laboratory, Pasadena, California, July 1962. 


where / is the frequency in hertz; v s is the shear velocity of 
the flow in meters per second, see [7]; x is the dimensionless 
frequency, x — 1,200 f/vm] and v m in meters per second 
is the mean velocity at 10 m. Denoting k = ( v s /v m ) 2 as 
the surface drag coefficient [Eq. (2)] gives 

«,(/) = 4, 800*. (jj-j^j) (3) 

The surface drag coefficient is obtained from the roughness 
of the terrain, see [7] 

k = (2.5/n(z/z 0 )) -2 (4) 

where z is the distance from the ground (z = 10 m) and 
is the roughness length (5-20 cm in the DSS-13 envi- 
ronment). For the roughness z Q — 0.05 m, one obtains 
k — 0.0057; and for = 0.2 m, k = 0.0105, which agrees 
with estimates by Levy and McGinness [7] of k at the an- 
tenna site in Goldstone, California. 

It is well known (see [8]) that the standard deviation of 
a stationary process can be determined from its spectrum, 
namely 

<7„ = (t„(oo)-T„( 0)) 1/2 (5a) 

where 

/ 

T v {f) = J S v (v)dv = 6Kv 2 m (l - (1 + s)- 1 ' 3 ) (5b) 

0 

thus, 

a v = av m , a = (6«) 1/2 (6) 

For a typical surface drag coefficient k — 0.005-0.010, one 
obtains a = 0.17-0.25. 

Although the force acts on the whole antenna struc- 
ture, in this article the wind force is assumed to act only 
on the antenna dish. This is a fair assumption, since most 
of the wind energy is absorbed by the tipping structure. 
The force distribution on the dish surface yields impor- 
tant information for determining the wind action on the 
antenna. By using the available quasi-static wind-tunnel 
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data, 6,7 the spatial force distribution can be obtained from 
the JPL-IDEAS model. 8 As a result, the wind force as a 
time and spatial variable becomes a time variable in the 
dynamic model. In the antenna model, the force acting 
on the dish has two components: F x and F y (wind from 
x and y directions). Since the time characteristics of the 
wind force do not depend on its direction, it is not nec- 
essary to distinguish the direction, thus the wind force is 
simply denoted by F. 

Statistical properties of the wind force are determined 
from the previously described wind- velocity parameters. 
At a given height, the force is proportional to the square 
of the velocity [9] 

F = kv 2 (7) 

and the velocity consists of a constant mean flow with 
a superimposed fluctuation, as in Eq. (1). Therefore, the 
force is decomposed into a steady force F m and a turbulent 
force F t with zero mean value 

F = F m + F t (8) 

The static force F m (lb) is determined for the speed v m 
(mph) 

Fm = F(v m ) = hv^ = aiF om (9) 

where a w = v m /100 , F om = 1 lb is the force at a wind 
velocity of 100 mph, and k — 0.0001F om . The gust F t , on 
the other hand, is determined as follows: 

F t = afVt (10) 

where aj = dF/dv \ v=Vm = 2 kv m = 2F m /v m depends 
on the wind velocity v m . From Eq. (10) the wind-force 
spectrum can be obtained as follows: 

Sj(f) = a)S v (f) (11) 

Thus, the standard deviation crj of the gust is 

<r f = \a f \<T v ( 12 ) 

6 R. B. Blaylock, op. cit. 

7 N. L. Fox, op. cit. 

8 R. Levy and D. Strain, op. cit 


Finally, by introducing Eq. (6) to Eq. (12), one obtains 

( T/ ~ 2 aF m (13) 

for the static force F m given by Eq. (9). 

III. Wind-Force Simulations 

Static and dynamic wind forces have been simulated. 
The static force is determined from Eq. (9), and the static 
force simulations are obvious. The dynamic force is gen- 
erated by a filter with a transfer function T(f ) 

T( n - a ° + + + ' ’ ’ + (ia\ 

b 0 + &is + &2$ 2 + • • • + b m s m 

The parameters a x and 6 X , i = 1, . . . ,m, m > 0 are de- 
termined so that the power spectrum of T(f) is the best 
approximation of the wind-power spectrum 5/(/). For 
this purpose, the square root sj of the power spectrum Sj 
is used 

S/if) =(s,(fj) 1/2 (15) 

The order m of the filter and the parameters a x and , i — 
l,...,m can be determined. For the frequency interval 
[/i, / 2 ], /2 > fi , denote the norm || . \\ w of a function g(f) 
as 



where w(f) > 0, fe [fi , / 2 ] is a weighting function. Denote 
the poles of the transfer function A,-, i = 1, . . . , m, then the 
filter parameters, and the filter order should be found so 
that the index J, defined as follows, 

J = || .,(/) - \T(f)\ ||. (17) 

is minimal, subject to constraints 

Re(\i) < 0, for i = 1, . . . ,m (18) 

The constraints assure filter stability. 
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The plot of the function sj(f ), obtained from Eq. (15) 
(the Davenport model for 30-mph wind speed), is shown 
by the solid line in Fig. 1. A linear filter that gener- 
ates the wind force of a spectrum close to the Daven- 
port model is determined by minimizing J subject to con- 
straints in Eq. (18). The values of the Davenport spectrum 
for high frequencies (in the case considered for frequencies 
1-1,000 Hz) are small; therefore, they should be weighted 
significantly to obtain a good fit in this frequency range. 
The following weights have been assigned: w = 3 for 
frequencies 0.1-1 Hz, w — 10 for 1-10 Hz, it; = 30 for 
10-100 Hz, and w = 100 for 100-1,000 Hz. For different 
integers m, the index J has been minimized; the integer 
m = 4 was found to be the lowest filter order for which 
the results have been satisfactory. The filter parameters 
for the minimal solution are: 

a 0 = 0.00037, a! = -0.92431, a 2 = -138.7468, 
a 3 = -56.4776, a 4 = -0.03952 
6 0 = —0.000306, bi = -0.18103, b 2 = -8.7769, 
b 3 = -57.9502, b 4 = -11.8569 

The magnitude of the filter- transfer function is shown 
by the dashed line in Fig. 1. The figure shows a good 
approximation of the Davenport spectrum within the in- 
terval of interest [fi = 0.001, / 2 = 10] Hz. The wind force 
generated by the filter for 30-mph wind is shown in Fig. 2. 

A similar approach was used to determine wind filters 
and wind forces for wind velocities v m = 10, 20, 40, and 
50 mph. 

IV. Pointing-Error Simulations 

The antenna state-space model with the position loop 
closed, described in [4], is as follows: 

x — Ax- b B c u + B W F , y = Cx 4- Du (19) 

where x is a state vector n x l;u is the command signal 
p x 1; F is the wind force (single input); y is the output 
q x 1. The input matrix B w is determined so that for 
a wind of 100 mph, a static wind force F om = 1 lb is 
applied. The plot of drag force versus antenna- elevation 
angle obtained from the wind-tunnel experiment 9 is shown 
in Fig. 3. In this article, two elevation positions of the 


9 Engineering Data 34- M Antenna, JPL Drawing 9493280 (internal 
document), Jet Propulsion Laboratory, Pasadena, California, 1988. 


antenna are simulated: 90 deg and 120 deg (that is, 60 deg 
with wind blowing from behind the dish). Note from Fig. 3 
that the drag is almost the same for the both positions. 
Nevertheless, it will be shown later that the effects of the 
wind actions on the antenna are distinct. 

Since the antenna model is linear, the static force (mean 
wind force) and dynamic force (gust force) are simulated 
separately. However, they can be superimposed, if nec- 
essary. For the static wind disturbances the system with 
the position loop closed [represented by Eq. (19)] is simu- 
lated with the zero command signal u and constant wind 
disturbance F m . The steady-state response (elevation and 
cross-elevation pointing error) is obtained for a 60- and 
90-deg antenna-elevation angle. The plots of pointing er- 
ror due to static wind from the x and y directions for differ- 
ent wind speeds are shown in Fig. 4. The following results 
for static wind load are obtained: For 30-mph x-direction 
wind, the elevation pointing error is 9.0 mdeg for a 60-deg 
elevation angle and 1.7 mdeg for a 90-deg elevation angle. 
The cross-elevation pointing error is 5.8 mdeg for a 60-deg 
elevation angle and 6.6 mdeg for a 90-deg elevation angle. 
For 30-mph y-direction wind, the elevation pointing error 
is 1.1 mdeg for a 60-deg elevation angle and 1.3 mdeg for 
a 90-deg elevation angle; the cross-elevation pointing error 
is 0.13 mdeg for a 60-deg elevation angle and 0.22 mdeg for 
a 90-deg elevation angle. The small pointing error due to 
wind in the y direction is a result of the integrating action 
of the controller. 

The gust wind load was applied to the antenna model. 
The responses of the antenna (elevation and cross- 
elevation pointing errors, elevation and azimuth encoder 
readings) are presented in Fig. 5 for x-direction wind and 
in Fig. 6 for y-direction wind. Wind force was simulated 
for velocities 10, 20, 30, 40, and 50 mph for 60- and 90-deg 
elevation angles of the antenna. The results are presented 
in Fig. 7 for x-direction wind and 60-deg elevation angle 
(solid line) and 90-deg elevation angle (dashed line). Fig- 
ure 8 shows the results for y-direction wind. Note that 
like the results of the static load, the results of the dy- 
namic loading show the pointing errors to be proportional 
to the square of the velocity (the symbols “o” in Figs. 7 
and 8 denote the simulation results, and the lines denote 
the least-squares quadratic approximation). It can be seen 
in Fig. 8 that the integrating action of the controller an- 
nihilated a significant part of the pointing error due to 
static y-direction wind; however, the dynamic part of the 
y-direction wind is not compensated for by the controller, 
which causes a significant pointing error. 

The spectra of the wind-gust response were analyzed. 
The spectra of the elevation and cross-elevation pointing 
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errors, and elevation and azimuth encoder readings for 
x-direction wind gusts are shown in Fig. 9, and in Fig. 10 
for y-direction wind gusts. The spectra show that modes of 
low frequency (2 Hz), medium (4 Hz), and high frequency 
(9 Hz) are important in wind-gust simulations. 

Finally, a power spectral-density analysis of the point- 
ing errors was performed for the antenna at the 90-deg 
elevation position, and results were compared with the 
simulation results of Figs. 5, 6, 9, and 10. The spectral 
analysis provides a precise measure of root-mean-square 
(rms) error, limited only by step-size effects. It serves as 
a convenient method for statistical verification of simula- 
tions. 

The power spectral-density function of the pointing er- 
ror S € (f) was calculated according to 

s c (f) = w 2 (f) S,(f) (20) 

where W(f) is the amplitude- versus- frequency response 
function between the respective wind-disturbance inputs 
and pointing-error outputs of the system described by 
Eq. (19), and S/(/) represents the wind-force spectrum 
from Eq. (11). The corresponding spectral distributions 
and rms errors were determined from a numerical integra- 
tion of S e (f) over frequency. An integration routine that 
adopted step size in the vicinity of peaks of S € (f) was 
employed to minimize computational errors. 


The frequency-response functions for elevation- and 
cross-elevation pointing errors are plotted in Figs. 11(a) 
and (b) for x-axis, and in Figs. 12(a) and (b) for y-axis 
wind direction. The corresponding power spectral density 
and spectral distribution functions for pointing errors are 
shown in Fig. 13 for x-axis winds and in Fig. 14 for y-axis 
winds. The mean wind velocity in both cases was 30 mph. 
The rms pointing errors were calculated from the corre- 
sponding spectral distribution functions at a frequency of 
10.0 Hz. For x-axis direction wind the elevation pointing 
error is 1.60 mdeg and the cross-elevation pointing error 
is 5.15 mdeg, while for y-axis direction wind, the elevation 
pointing error is 5.53 mdeg, and the cross-elevation point- 
ing error is 0.65 mdeg. These errors agree well with the 
values obtained from simulations of Figs. 5 and 6. 


V. Conclusions 

Wind forces acting on the antenna structure have been 
modeled in this article. Their static and dynamic parts 
have been applied as a disturbance to the closed-loop an- 
tenna model at 60- and 90-deg elevation angles. The re- 
sults show that both the static and dynamic parts of the 
pointing error are of equal importance, and that despite 
the dominant presence of low-frequency components in the 
wind-disturbance spectrum, the high-frequency modes of 
the antenna have been excited, which caused significant 
pointing errors. In order to improve the analytical results, 
experimental verification should be the first priority. 
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Fig. 1. Square root of power spectrum of the Davenport wind-gust model, for 30-mph wind at 

the DSS-13 antenna site. 



Fig. 2. Wind force generated by the filter for a 30-mph wind gust. 
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Fig. 4. Pointing error due to static wind from x - and y-direction for different wind speeds: 
(a) elevation pointing error; (b) cross-elevation pointing error; (c) elevation encoder reading; 
and (d) azimuth encoder reading. 
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Fig. 9 (contd). 
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Fig. 10 (contd). 
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